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Abstract 

I review the problem of parton saturation and its implications through three in 
principal different approaches, but somewhat related: saturation in a geometrical 
approach, QCD saturation through the Color Glass Condensate and perturbative 
Pomeron approach with initial conditions. I also make some comments about how 
that could be related to a Quark Gluon Plasma formation. 
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1 Introduction 

In the recent experiments like DIS at HERA or the heavy-ion experiments at RHIC, 
and also in expected LHC at CERN, the number of involved partons is very large, 
due to the high energy and/or the high number of participants of those experiments. 
These high parton densities should in principal lead to an extremely huge multiparticle 
production, but experimentally we have seen that this is not the case. So there should 
be a mechanism that reduces the number of created particle. Two kind of phenomena 
have been proposed. On one side there is the possibility of saturation in the initial 
state of the collision. On the othe side, there is the proposal of the creation of a Quark 
Gluon Plasma, and it has been presented as a final state interaction mechanism. Here, 

1 review the problem of parton saturation and its implications through three in princi- 
pal different approaches, but somewhat related: saturation in a geometrical approach, 
QCD saturation through the Color Glass Condensate and perturbative Pomeron ap- 
proach with initial conditions. I will also make some comments about how that could 
be related to a QGP formation. 

2 Geometrical approach to saturation 
2.1 String models and percolation 

In many models of hadronic and nuclear collisions, color strings are exchanged between 
the projectile and the target. Those strings act as color sources of particles through 
the creation of q — q pairs from the sea. The number of strings grows with the energy 
and with the number of nucleons of the participant nuclei. 

In impact parameter space these strings are seen as circles inside the total collision 
area. When the density of strings becomes high the string color fields begin to overlap 
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and eventually individual strings may fuse, forming a new string -cluster- which has 
a higher color charge at its ends, corresponding to the summation of the color charges 
located at the ends of the original strings. The new string clusters break into hadrons 
according to their higher color. As a result, there is a reduction of the total multiplicity. 
Also, as the energy-momenta of the original strings are summed to obtain the energy- 
momentum of the resulting cluster, the mean transverse momentum of the particles 
created by those clusters is increased compared to the one of the particles created from 
individual sources. 

As the number of strings increases, more strings overlap. Some years ago, it has 
been proposed in Ref. pQ that above a critical density of strings percolation occurs, so 
that paths of overlapping circles are formed through the whole collision it is 

represented in Fig. [TJ Along these paths the medium behaves like a color conductor. 
Also in pQ, we have made the remark that several fused strings can be considered as 
a domain of a non thermalized Quark Gluon Plasma. The percolation gives rise to the 
formation of a non thermalized Quark Gluon Plasma on a nuclear scale. 

Note that here we are not speaking about a final state interaction phenomenon, since 
there is no thermalization involved. In fact, what we are trying is to determine under 
what conditions the initial state configurations can lead to color connection, and more 
specifically, if variations of the initial state can lead to a transition from disconnected 
to connected color clusters. The results of such a study of the pre-equilibrium state 
in nuclear collisions do not depend on the subsequent evolution and thus in particular 
not require any kind of thermalization. 

The structural problem underlying the transition from disconnected to connected 
systems of many components is a very general one, ranging from clustering in spin 
systems to the formation of galaxies. The formalism is given by percolation theory, 
which describes geometric critical behavior. 
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2.2 Percolation theory 

Consider placing N small circular discs (color sources, strings or partons) of radius 
r onto a large circular manifold (the transverse nuclear plane) of radius R; the small 
discs may overlap. With increasing density 

Nirr 2 

this overlap will lead to more and larger connected clusters. The most interesting 
feature of this phenomenon is that the average cluster size increases very suddenly 
from very small to very large values. This suggests some kind of geometric critical 
behavior. In fact, the cluster size diverges at a critical threshold value r] c of the density. 
This appearance of an infinite cluster at rj = r] c is defined as percolation: the size of 
the cluster reaches the size of the system. rj c has been computed using Monte Carlo 
simulation, direct connectedness expansion and other different methods. All the results 
are in the range 7] c = 1.12 -f- 1.175. 

In our model [l] we had proposed a fixed radius for the independent color sources 
of r = 0.2 -r- 0.25 fm, that corresponds to a momentum around 1 GeV. This value has 
been obtained from Monte Carlo simulations in the framework of the String Fusion 
Model Code (SFMC) [2] made at SPS energies. According to eq. (JTJ), in order to 
estimate the density rj, one needs to know the number of sources N. In our model, it 
is obtained from the SFMC, that, for nucleus-nucleus collisions, takes into account two 
contributions: one proportional to the number of participant nucleons -valence-like 
contribution- and another one proportional to the number of inelastic nucleon-nucleon 
collisions. Note that N will depend on the energy y/s (or equivalently, on x) and on the 
number of participant nucleons A, so in some way the condition to achieve percolation 
depends on A and s, rj = rj(A, x). ttR 2 corresponds simply to the nuclear overlap area, 
Sa, at the given impact parameter. This overlap area can be determined in a Glauber 



study, using Woods-Saxon nuclear profiles. That leads to the following results: In our 
model, at SPS energies, the critical threshold for percolation could eventually been 
achieved for the most central Pb-Pb collisions, and for sure in Au-Au central collisions 
at RHIC energies and even in p-p collisions at LHC energies. Just as an example, 
in first approximation, one can estimates analytically that at very high energies, for 
central A-A collisions, if we take the number of initially created sources as proportional 
to the number of collisions, N cx v4 4//3 and the nuclear overlap area Sa = nR 2 oc A 2 ^ 3 
then 

_ Nirr 2 _ N 8aurce8 Si ^ a2/3 ^ 

ttR 2 Sa 

2.3 The size of the color sources 

Let us now introduce a new question: If you look at a fast nucleon coming at you, what 
do you see? It depends on who is looking. Another nucleon sees a disc of radius r ~ 1 
fm and a certain greyness. A hard photon, with a resolution scale Q^ 1 << 1 fm, sees 
a swarm of partons. How many there are depends on the resolution scale: given a finer 
scale, you can see smaller partons, and there are more the harder you look. The partons 
in a nucleon have a transverse size determined by their transverse momentum kx , 
with Tt ~ 1/kx ■ The scale Q^ 1 specifies the minimum k^ 1 resolved, so the probing 
photon sees all partons in the range < kx < Q, or equivalently it sees all the partons 
with a radius tt > Q -1 ■ 

So the partonic size is through the uncertainty relation determined by its average 
transverse momentum, r 2 ~ 1/ < k^ >, for a given resolution scale, < k\ >~ Q 2 . 

These ideas are illustrated in Figs. El and El 

In order to know if the percolation density is achieved in a collision, we need to 
compute our eq. (JTJ). Then it is necessary to know the number of initially created 
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sources (partons or strings, it depends on the kind of model we are using) N, and the 
size of those sources r. Remember that in the previous section we have used a fixed 
size of r = 0.2 4- 0.25 fm that corresponds to a moment of 1 GeV. In pQ N is calculated 
in the framework of the SFMC. 

In the number of sources N is calculated through the Wonded Nucleon Model 
for nucleus-nucleus collisions, using P.D.F.'s for the nucleonic density. The authors use 
here an slightly different strategy: Instead of fixing the radius of the initially created 
sources, they estimate the momentum of those sources that will lead to percolation. 
They transform then eq. (JIJ into the following one: 

iWr 2 N 

Vc = 7# = gpp ' (3) 

The condition for percolation, taking rj c = 1.12 will be then: 

Q 2 = — (4) 

They find that for central Pb-Pb collisions at SPS energies, the critical momentum for 
percolation is Q 2 C ~ 1 GeV 2 , in accordance with our result. For central Au-Au collisions 
at RHIC, Ql ~ 2.5 GeV 2 . Note that again the condition to achieve percolation depends 
on A and y/s, Q c = Q c (A,x). 

Beyond the percolation point, one has a condensate, containing interacting and 
hence color-connected sources of all scales k T < Q. The percolation point thus specifies 
the onset of color deconfinement; it says nothing about any subsequent thermalization. 



3 The Color Glass Condensate 

Now we arrive to another approach, the QCD saturation through the formation of a 
Color Glass Condensate (CGC) 0]. The idea is the following: At high energy, the QCD 
cross-sections are controlled by small longitudinal momentum gluons in the hadron 
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wave function, whose density grows rapidly with increasing energy or decreasing x, 
due to the enhancement of radiative process. If one applies perturbation theory to this 
regime, one finds that, by resumming dominant radiative corrections at high energy, 
the BFKL equation leads to a gluon density that grows like a power of s and in 
consequence to a cross-section that violates the Froissart bound. Nevertheless, the use 
of perturbation theory to high-energy problems is not obvious. In fact, the BFKL and 
DGLAP equations are linear equations that neglet the interaction among the gluons. 
With increasing energy, recombination effects -that are non-linear- favored by the high 
density of gluons should become more important and lead to an eventual saturation of 
parton densities. 

These effects become important when the interaction probability for the gluons 
becomes of order one. Taking as the transverse size of the gluon and xG< ^^ - as 
the density of gluons, the interaction probability is expressed by 



Equivalently, for a given energy, saturation occurs for those gluons having a sufficiently 
large transverse size r\ ~ l/Q 2 , larger than some critical value l/Q s (x,A). So the 
phenomenon of saturation introduces a characteristic momentum scale, the saturation 
momentum Q s (x,A), which is a measure of the density of the saturated gluons, and 
grows rapidly with 1/x and A (the atomic number). The probability of interaction - 
that can be understood as "overlapping" of the gluons in the transverse space- becomes 
of order one for those gluons with momenta Q 2 ^ Q s (x, A) where 



For Q 2 < Q 2 (x,A), the non-linear effects are essential, since they are expected 
to soften the growth of the gluon distribution with r = ln(l/x). For a nucleus, 
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a s N c 

Q 2 



X 



xG(x,Q 2 ) 
^R 2 



(5) 





xGa(x, Ql) oc A and irR A oc A 2 / 3 , so eq. (0) predicts Q 2 S oc A 1 / 3 . One can esti- 
mate the saturation scale by inserting the BFKL approximation into eq. (jHJ). This 
gives (with 5 ~ 1/3 and A ~ ca^ in a first approximation): 

Qfc.A) ~ A 5 x~\ (7) 

which indicates that an efficient way to create a high-density environment is to combine 
large nuclei with moderately small values of x, as it is done at RHIC. In fact the 
estimated momentum for saturation at RHIC will be Q s = 1 -j- 2 GeV, in accordance 
with the result of the previous section. 

This equation also shows that for sufficiently large energy or x small enough, 
Ql(x,A) ^> Aq CD and a s (Q s ) <C 1, which characterizes the regime of weak coupling 
QCD. But although the coupling is small, the effects of the interactions are amplified 
by the large gluon density: at saturation, Ga(x,Q^) ~ l/a s (Q s ) ^> 1, so the gluon 
modes have large occupation numbers, of order l/a s (corresponding to strong classical 
fields A ~ l/g), which suggests the use of semi-classical methods. 

Note that in the string models described above, the formation of new strings - 
sources- through the creation of q — q pairs from the sea is a non-perturbative way of 
gluon emission. 

One can write a classical effective theory based on this general idea: the fast par- 
tons -valence quarks with large longitudinal momentum- are considered as a classical 
source that emits soft gluons -with smaller longitudinal momenta- which are treated 
as classical color fields A[p]. The fast partons move nearly at the speed of light in the 
positive x + direction, and generate a color current = <5 M+ p. By Lorentz contraction, 
the support of the charge density p is concentrated near the light-cone longitudinal 
coordinate x~ = 0. By Lorentz time dilation, p is independent of the light-cone time 
x+. 
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The Yang Mills equation describing the soft gluon dynamics reads 

Dv p»v = s» + p(x-,x) . (8) 

Physical quantities, as the unintegrated gluon distribution, are obtained as an average 
over p: 

< A\X) A\Y) > X =J D[p] W x [p] A*[p}(X) A l [p]{Y) , (9) 

where A l (X) corresponds to the classical solution for a given p, and W x {p] is a gauge- 
invariant weight function for p. What we are doing is a kind of Born-Oppenhaimer 
approximation: first, we study the dynamics of the classical fields (Weizsacher- William 
fields) for a given configuration p of the color charges, and second, we average over all 
possible configurations. 

For the classical solution we find: 

F +l (x~,x ± ) = 6(x-)-V(x ± )(d i V(x ± y) = d + A i , (10) 
9 

where V(x±) is the Wilson line 

V ] {x±) = Pexp{igJdx'A + {x-,x)} (11) 

and A + [p] is the solution of the equation of motion (jHJ) in the covariant gauge: —V\A + = 
P- 

The weight function W T [p] is obtained by integrating out the fast partons, so it 
depends upon the rapidity scale r = ln(l/a;) at which one considers the effective 
theory This can be taken into account via a one-loop background field calculation, 
and leads to a renormalization group equation (REG) for W T [p] which shows how the 
correlations of p change with increasing r. Schematically: 

dW T [p] 



Or 



\LA< [A m wm - (12) 



This is a functional diffusion equation, where the kernel x[p\ plays the role of the dif- 
fusion coefficient in the functional space spanned by p(x~,x±). This equation encom- 
passes previous evolution equations [5] developed by Balitsky, Korchegov and Weigert. 
This kernel is positive definite and non-linear in p to all orders. It depends upon p via 
the Wilson line (fTT|). 

As it has been said above, the physical quantities, as the gluon density, are obtained 
as an average over p: 

n{x, k ± ) = -L {F +\k ± )F+*(-k ± )) x . (13) 

In order to calculate this average, we need to solve the REG (|12p. Approximate 
solutions to this equation can be obtained in two limiting cases: 

• At low energy, or large transverse momenta k\ ^> Q 2 s (x), we are in a dilute regime 
where fields and sources are weak, and the Wilson lines can be expanded to lowest order 
in A + , V'(x±) ~ 1 + ig j dx~A + (x~, x±). In this case, the REG equation reduces to 
the BFKL equation, and the gluon density of eq. (fT3*|) grows both with l/k\ and 1/x 
(Bremsstrahlung): n(x,k±) ~ ^- 

• At high energies, or low momenta k\ ^ Q 2 (x), the color fields are strong, A + ~ 
1/g, so the Wilson lines rapidly oscillate and average away to zero: V ~ ~ 0. Then 
the kernel x becomes independent of p, and we obtained a gluon density that increases 
linearly with the evolution "time" r = ln(l/x): n(x,k±) ~ _L \ n ^ j n I_ That 
is, we find saturation for the gluon density, that grows logaritmically with the energy 
since r ~ In s: unitarity is restored. 

We call the high density gluonic matter at small- a; described by this effective theory 
a Color Glass Condensate: Color since gluons carry color under SU(N C ); Glass since 
we have classical coherent fields which are frozen over the typical time scales for high- 
energy scattering, but randomly changing over larger time scales. So we have a random 
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distribution of time-independent color charges which is averaged over in the calculation 
of physical observables, in order to have a gauge independent formulation -in analogy 
to spin glasses-; and Condensate because at saturation the gluon density is of order 
l/a s , typical of condensates, so we have a system of saturate gluons that is a Bose 
condensate. 

4 Perturbative QCD pomeron with saturation in 
the initial conditions 

Let us now try a different approach [6 . Consider now the nucleus-nucleus interaction 
as governed by the exchange of pomerons, as it is represented in Fig. |3J Its propagation 
is governed by the BFKL equation [7] . Its interaction is realized by the triple pomeron 
vertex. Equations which describe nucleus-nucleus interaction in the perturbative QCD 
framework have been obtained in [Hj. They are quite complicated and difficult to solve, 
but they will be not needed for our purpose. Knowing that the AGK rules are satisfied 
for the diagrams with BFKL pomerons interacting via the triple pomeron vertex [Oj 
and using arguments of [Ej it is easy to conclude that the inclusive cross-section will 
be given by the convolution of two sums of fan diagrams propagating from the emitted 
particle towards the two nuclei. 

Taking A = B and constant nuclear density for |6| < R A , one can find the inclusive 
cross-section in perturbative QCD as 

I A (y,k) = A^TtRf-^ J d 2 re* kr [A$ A (Y-y,r)][A$ A (y,r)], (14) 

where A is the two-dimensional Laplacian and $(y, r) is the sum of all fan diagrams con- 
necting the pomeron at rapidity y and of the transverse dimension r with the colliding 
nuclei, one at rest and the other at rapidity Y. The function <p A (y, r) = $(?/, r)/(27rr 2 ), 
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in the momentum space, satisfies the well-known non-linear BK equation [5J 

^ = -%?)-^ l g) 1 (is) 

dy 

where y = dy, a = a s N c /n, a s and N c are the strong coupling constant and the 
number of colors, respectively, and H is the BFKL Hamiltonian. This equation has to 
be solved with the initial condition at y = determined by the color dipole distribution 
in the nucleon smeared by the profile function of the nucleus. 

We can take the initial condition in accordance with the Golec-Biernat distribution 
[12], that takes into account saturation according to the CGC, duly generalized for the 
nucleus: 

0(0, q ) = --a Ei ( - ? ) , (16) 

rv ,qj 2 I 0.3567 GeV 2 / ' V ; 



with 



a = A*™%*. (17) 



Then we can find the following result: One can observe that whereas at relatively 
small momenta the inclusive cross-sections are proportional to A, that is to the number 
of participants, at larger momenta they grow with A faster, however noticeably slowlier 
than the number of collisions, approximately as A 1A . The interval of momenta for which 
I a oc A is growing with energy, so that one may conjecture that at infinite energies all 
the spectrum will be proportional to A. 

One can also compare ^3] the results obtained on the basis of the expression for 
the inclusive cross-section which follows from the AGK rules applied to the diagrams 
with QCD pomerons interacting via the three-pomeron coupling [8 with a somewhat 
different expression obtained from the color dipole picture [TJ]. These two formally 
different expressions lead to identical results at momenta of the order or higher than 
the value of the so-called saturation momentum Q s . At momenta substantially lower 
than Q s the color dipole cross-sections differ from the ones from the AGK rules by a 
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universal constant factor < 0.8 4- 0.9. In both cases the spectra at momenta below Q s 
are found to be proportional to the number of participants (A for collisions of identical 
nuclei) and not to the number of collisions v4 4//3 . Since Q s grows with energy very fast, 
the region where the spectra are proportional to A extends with energy to include all 
momenta of interest. At momenta greater than Q s the spectra grow faster than A but 
still much slowlier than v4 4//3 (a numerical fit gives something like A 1 - 1 ). 

Note that in the string models, like Quark Gluon String Model or Dual Parton 
Model ^2] each string corresponds to the exchange of pomerons, and the interaction 
among the strings would correspond to interaction among the pomerons through the 
triple pomeron vertex. In particular, in the DPM some kind of saturation is also 
included through the shadowing in the initial state of the collision JH]. in fact, in 
absence of shadowing, the DPM at high energies leads to multiplicities that scale with 
the number of binary collisions rather than to a scaling with the number of participants. 
This is a general property of Gribov's Reggeon Field Theory which is known as AGK 
cancellation, analogous to the factorization theorem in perturbative QCD and valid for 
soft collisions in the absence of triple Pomeron diagrams. Because of this, a dynamical, 
non-linear shadowing has been included in the DPM. It is determined in terms of the 
diffractive cross-section. It is controlled by triple pomeron diagrams, shown in Fig. El 
and it should lead to saturation as s — > oo. In fact, the shadowing corrections in this 
approach have a positive contribution to diffraction and a negative one to the total 
cross-section. 

5 Conclusions 

We have compared different models that takes into account saturation in different 
ways: from the semi-phenomenological fusing color sources picture for the soft domain 
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including percolation, the QCD saturation through the Color Glass Condensate and 
those which follow from the pomeron approach, perturbatively derived from QCD, 
taking into account saturation in the initial conditions. 

In fact, it seems that the exchanged of elemental objects, color sources -strings, 
partons or pomerons-, should lead to a saturation in the initial conditions when the 
densities are high enough. 

I would like to finish making some remarks. Concerning the relation between per- 
colation and thermal phase transitions, some thermal critical behavior, such as the 
magnetization transition for ferromagnetic spin systems, can be equivalently formu- 
lated as percolation. However, percolation seems to be a more general phenomenon 
and in particular can occur even when there is no thermal critical behavior. A spe- 
cific example of this is the Ising model in a non-vanishing external field, which has a 
percolation transition even though there is no magnetization transition. 

In nuclear collisions there is indeed, as function of parton or string density, a sudden 
onset of large-scale color connection. There is a critical density at which the elemental 
objects form one large cluster, losing their independent existence and their relation to 
the parent nucleons. Percolation would correspond to the onset of color deconfinement 
and although it may be a prerequisite for any subsequent QGP formation, it does not 
require or imply any kind of thermalization. 

The ancient idea about the formation of an ideal QGP -weakly coupled plasma of 
quarks and gluons- above a critical temperature T c ~ 160 MeV has changed: It seems 
that for moderate temperatures (1 -i-3) T c , accessible at RHIC, the plasma is predicted 
by non-perturbative lattice to be strongly coupled (sQCD), and that its source could 
be a saturated initial state [Ej of the type I have described here. 
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Figure 1: From left to righ: Density of strings in the transverse space, from low energy 
and/or low number of participants to high energies and/or high number of participants. 
In the last circle we show percolation. 
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Figure 2: From left to righ: The structure of an incoming nucleon seen for increasing 
resolution. 



18 




19 



Figure 4: A typical diagram for the inclusive cross-section in nucleus-nucleus collisions. 
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Figure 5: Fan diagram with triple pomeron interaction. 
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